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1 Introduction 

Nuclear shadowing or the depletion at low Bjorken x of the nuclear Deep Inelas- 
tic (DI) structure function, F^, with respect to the nucleon one, , has been 
observed in a number of experiments (seeEJ and references therein). Assuming 
the universality of parton distributions in nuclei, one expects nuclear shadow- 
ing to be present in other high energy processes as well, such as Drell-Yan pair, 
J/jp and T production in lepton- nucleus, hadron-nucleus and nucleus-nucleus 
collisions. In particular nuclear shadowing might be concurring with the other 
mechanisms among which quark-gluon plasma formation, rthat result in a de- 
pletion of the observed cross sections for these processes Ha. A quantitative 
understanding of both the x and dependences of the nuclear parton distri- 
butions at low and very low x is therefore a necessary step for interpreting the 
outcome of future experiments at RHIC and at the LHC. Recent calculations 
rely on non-perturbative models for the nuclear parton distributions at a given 
(low) scale, Q^, combined with DGLAPcl perturbative evolution. They are all 
therefore affected by the uncertainty in the initial parton distributions and, 
in particular, in the gluon distribution which governs evolution at low x, and 
which is poorly known experimentally. Moreover, evolution to moderate values 
of {Q^ ~ lOGeV^) depends even more strongly, by the value of Ql itself 
which, if varied within the range, Ql = 0.8 — few GeV^ can lead to completely 
different evolution patterns for the shadowing of both the structure functions 
and the gluon distribution. These simple facts hamper the possibility of pre- 
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dieting the dependenee of nuelear shadowing although its behavior within 
pQCD is quite weU understood and it can be predicted accurately for each 
choice of initial conditions. 

In this paper we suggest an avenue for dealing with the model depen- 
dence associated with the initial conditions of perturbative evolution, which 
is based on a study of the low x and high asymptotic behavior of the 
shadowing ratios Rq = Ga/Gn and Rp = F^/F^, Gjv(a) and F^^^^ be- 
ing the gluon distributions and the Deep Inelastic (DI) structure function in 
a nucleon, (N), and in an isoscalar nucleus, {A), respectively. In the asymp- 
totic regime defined by and x ^ Xo, Xq < 0.1, tbe Double Log- 
arithmic Approximation (DLA) to pQCD evolution jagplies The proton 
structure function data analyzed recently at HERA eIlI have been proven U 
after the suggestion of Refl3, to evolve according to DLA. The key test is 
to show that the data obey Double Asymtptotic Scaling (DAS) in the vari- 
ables p = 7((r - Yo)/0^/^ and a = j-^iY - Yo)^^/^, 7 = 6/(33 - 2Nf)^/^, 
Y = Inl/x, C = "/^HIuQ^/AIcd/'^^QI/^qcd)- Violations from DAS (other 
than due to the fact, that the data lie in a pre-asymptotic region or to NLO 
DGLAP corrections lab) would signal either the onset of other approximations 
such as leading- (next-to-leading)-log in 1/x, LLx (NLx) resummation, or the 
beginning of parton saturation. Here we wish to apply a DAS type of analysis 
to nuclear DI scattering. 

As a first step we show that if evolution proceeds through ordinary asymp- 
totic pQCD evolution equations in nuclei as well as in a free nucleon, the ratios 
Rp and Rq become a function of p only, which is per se a model independent 
result. We use this result as a basis for exploring the origin of scaling viola- 
tions in nuclei. In nuclei in fact, the asymptotic regime is entered in principle 
at different values of x^ and Q^^than in the proton. In particular Unitar- 
ity Shadowing Corrections (USC) EHlj'EiEj are expected to affect evolution at 
x^ > Xo because of the increase of the gluon density in a nucleus due to the 
overlapping of nucleons in the longitudinal direction. On a more speculative 
basis one might also expect the transition to the ln(l/a;) resummation to ap- 
pear in a different regime. In our approach such questions can be addressed 
systematically as they introduce specific scaling violations from the DLA result, 
appearing as a cr dependence in the ratios Rq and Rp. 

In summary, although it is technically predictable that in a proton at 
very low values of x and sufficiently large Q^, i.e. deeply in the asymptotic 
region, DGLAP evolution and the DLA should break down and give way to 
ln(l/a;) summation and to USC, it is still a major task to be able to pinpoint 
where and if the transition from the different regimes is going to take place 
in the kinematical regimes currently under exploration. Our goal is to obtain 
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some new insight by using nuclear targets where the asymptotic regime can 
in principle be reached at larger x. As a by-product we obtain quantitative 
predictions for RHIC and the LHC. 



2 DGLAP Evolution in Nuclei 

We first summarize results for ordinary DGLAP evolution applied to the nu- 
clear ratios at low x, assuming that the proton and the nuclear distributions 
evolve along similar paths. As it is well known evolution is driven by the gluon 
distribution which dominates over the sea quarks one and one can predict the 
behavior of the shadowing ratios, Rq and Rp with Q^: 

dRc „ f\ fx ,^2AG^v(y^ 

X [RG{y,Q')-RG{x,Q')] ^ 

y 

dRp f\ fx ,^2.\GNiy,Q^) 



f PgG(-,as{Q^) 

J X W 



X [RG{y,Q^)~RF{x,Q^)\ ^ , (2) 

y 

where we have disregarded the sea quarks distribution on the r.h.s. of the 
coupled DGLAP evolution equations; P^g and Pgg are the splitting func- 
tions evaluated at NLO; and we used the approximation F^^^^ w 5/181] jv(A)- 
Eqs.(l) and (2) show that the dependence of Rg and Rp is regulated by 
ia_subtle balance involving the parton distributions and the ratios themselves 
The following predictions can be made for the ratios Rg and Rp: If, 
as predicted by current non-perturbative shadowing models Rg is a growing 
function of x, then it also grows with Q^, the r.h.s. of Eq.(l) being positive 
{y > x). On the other side, defining Rg{x,QI) = Rq and Rf{x,QI) = Rp, 
one obtains two opposite behaviors for Rp namely: if Rp < Rq, then Rp 
grows with Q^; ii) if instead Rp > Rq, then Rp initially decreases with 
until it reaches the value of R^ and it subsequently starts increasing along with 
Rg- □ The "rate" of change with is governed by the ratios GN{y)/GN{x) 
and GN{y)/^N{x), at = Ql, respectively. If Ql < IGeV^, then a rapid 

"^Although the form of Eq.(l) might seem suggestive of a fixed point behavior Q, this is a 
priori not the case, since the quantity dGA{x, Q^)/d\nQ-^ /dGisfi^, Q^)/dlnQ^ depends on 
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evolution strongly reduces the shadowing in both Rq and Rp, between and 
« 2 — 3 GeV^ . If on the contrary ranges from 2 to 5 GeV^ then the 
evolution is slower and shadowing remains large even at « 10 — 100 GeV^ 
(these results affect quarkonia production where = Mj^^ w 10 GeV^ and 
Q2 EE Af2 w lOOGeV^). 



3 Double Asymptotic Scaling 

We now examine the nuclear DI structure function and gluon distribution 
in the asymptotic regime. The derivation of the equations of the DLA in a 
nucleus parallels the one for the proton, namely one first writes the LO DGLAP 
evolution equation for the gluon distribution in the limit x — > 0, in moments 
space (we have omitted the subscripts N{A) unless necessary): 

dgjn, Q^) as{Q^) o ^ ^ / (o^ 
-d]^ = ^^^OGin)9in,Q), (3) 

7gg(n) « 2Ca/{tT' ~ 1) being the anomalous dimension in the limit n ^ \. 
Solutions in (x, Q^) are found by evaluating the anti-Mellin transform, 

G(x,Q2) = / dng{n,Ql) 

J (J 

exp{(n-l)y + e/(n-l)}, (4) 

with the saddle point method. In Eq.(^ ff("', Qo) = x'^^''^^^ g{x,QQ), 

g{x, Q^) being the inital gluon distribution, and G{x, Q^) = xg{x, Q^). If one 
takes a "soft" initial condition such as, G{x,Qq) « A]\fX~^, A « 0, and Y 
and ^ are both similarly large, then g(r7-o,Qo) ~ ^N/in — 1), and the saddle 
point is: no - l/(2Ar) + ^l/(4Ar)2 + ^/AF w y^T+JjAY, AY ^Y- Y^, 
Yo = ln(l/a;o), yielding: 



with 'g{n,Q^) = g(n, Q^) exp[(n — 1)^]- ^ In terms of the DAS variables, p 
and a, no — I + pjy^ and: 

G = G^^-4(x,Q2)^/(p,^)exp(27fT), (6) 

where /(p, cr) is a function that depends on the explicit form of the initial 
conditions. The main result is that \n{G^^^)/ f{p, a)) is predicted to be linear 



We have omitted sub-leading corrections for simplicity. 
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in cr, the slope being fixed by 7. A similar behavior is found for F2 by solving 
the equation: 

dFi'{x,Q'^) 2a.c 



9 TT 



(7) 



SlnQ2 

Deviations fi'om DAS due to hard initial conditions i.e. when A > 0.2 are 
technically obtained as follows: g{no, Qq) « AN/{n — (A + 1)) and a physical 
solution in the limit F ^ is 

corresponding to the saddle point value, 710(A) — l/AY + A + 1 (see also&Q). 
Physically this behavior supports the appearance of USC which becornc more 
important just because of the steep rise at small x. The data of BpM favor 
the DAS scenario, whereas ±he more recent data at lower x and Q^U are best 
interpreted in terms of USCEEl. In this context it is therefore important to study 
nuclear shadowing because, due to the A-dependent increase in gluon density, 
the physics we are interested in is expected to show up at larger x values and 
because one can tune the variable A in order to discriminate among models. 

In a nucleus we make the two following assumptions for the asymptotic 
regime: i) the initial distributions are shadowed because of some non-perturbative 
mechanism, the evolution equations are not affected by the medium; ii) be- 
sides the non-perturbative shadowing the evolution equations have screening 
corrections at a larger xa. The derivation of the ratio Rq in the hypothesis i) 
yields: 
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where /a(p) = GA{no,Ql)/GN{no,Ql), and the functions /ij^^^^ are listed 
explicitely elsewhere. The main point is that the exponential terms appearing 
in G^^^, Eq.(H) calculated for a nucleus and for a nucleon respectively, cancel 
exactly and the remaining dependence on a is small asymptotically. This 
scaling result (shown in Fig.l) will be our reference point. 

The onset of a different evolution mechanism in the nucleus will appear as 
a a-scaling violation modifying the exponential behavior of Eq.(^). In general 
one can write 

Rg = Rg'"\p) X exp {2^{aA - a)} , (10) 
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Figure 1: Nuclear structure funcifeisins ratios, Rp, vs. the DAS variables p (top) and a 
(bottom). The theoretical curves E3 show the onset efi scaling in a. The non-perp™+>iitive 
shadowing pwdel used in the calculation is from Reflij. Experimental data fromtZHij (Li, 
C, Ca) and La (Xe). 



where i?^*^™ is the large a Umit of Eq.(^. The form of a a = aA{x,Q^) 
describes the approximation to pQCD evolution in a nucleus. Iruthe case of a 
"hard pomeron" boundary conditions plus screening correctionsEI it is 

aA=DAix,Q^)(^Xap+^^y (11) 

Da being a damping correction (more details are going to be found inE^I). 

In conclusion, we have shown model independent predictions for gluon 
shadowing in nuclei in the asymptotic region. We have also outlined an ap- 
proach for more detailed studies of the violations from DAS in nuclei as a 
signal of pQCD evolution mechanisms other than DGLAP or the DLA. Our 
calculations will be relevant in the regime accessible at, future experiments at 
RHIC, LHC and possibly at the eA project at DESYEII. 



Acknowledgments 

The authors wish to thank the Institute of Nuclear and Particle Physics at the 
University of Virginia for financial support. 

References 

[I] M. Arneodo, Phys. Rep. 240 (1994) 301. 

[2] Zi-wei Lin and M. Gyulassy, Phys.Rev.Lett. 77 (1996) 1222. 

[3] S. Liuti and R. Vogt, Phys. Rev. C51 (1995) 2244. 

[4] L.V. Gribov and Lipatov, Yad.Fiz. 20 (1975) 181; G. AltareUi and G. 

Parisi, Nucl.Phys. B126 (1977) 298; Dokshitzer, Sov.Phys.JETP 46 

(1977) 641. 

[5] A. De Rujula et al., Phys. Rev. DIO (1974) 1649. 

[6] HI Collaboration, T.Ahmed et al., Nucl.Phys. B439 (1995) 471; C. 
Adloff et al., Nucl.Phys. B497 (1997) 3; S. Aid et al., Nucl.Phys. B470 
(1996), 3. 

[7] ZEUS Collaboration, M. Derrick et al., Z.Phys. C69 (1996) 607; J. 

Breitweg et al., Phys.Lett. B407 (1997) 432 
[8] R. Ball and S. Forte, Phys. Lett. B335 (1994) 77; ibid B336 (1994) 77; 

Nucl. Phys. Proc. Suppl. 54A (1997) 163. 
[9] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1. 
[10] A.H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. (1994) 44. 

[II] J.C. Collins and I. Kwiecinski, Nucl. Phys. B335 (1990) 89. 
[12] K. Eskola, J. Qiu and X. Wang, Phys.Rev.Lett. 72 (1994) 36. 
[13] J. Qiu, Nucl. Phys. B291 (1987) 746. 

[14] Z. Huang, H.J. Lu and I. Sarcevic, Nucl. Phys. A637 (1998) 79. 

[15] L.L. Frankfurt and M.I. Strikman, Nucl. Phys. B316 (1989) 340. 

[16] E. Gotsman, E Levin, U.Maor and E. Naftah, Nucl. Phys. B (1999), 

|hep-ph/9808257 . 
[17] M. Arneodo et al., Nucl. Phys. B441 (1995) 12. 
[18] P. Amaudruz et al., Z. Phys. C53 (1992) 73. 

[19] M.R. Adams et al., Phys. Lett. B287 (1992) 375; Phys. Rev. Lett. 68 

(1992) 3266. 
[20] F. Cano and S. Liuti, in preparation. 

[21] Workshop on Physics with HERA as eA Collider, DESY, Hamburg, 25- 
26 May, 1999, C99/05/95. 



7 



